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Abstract 

A significant mathematical error is identified and corrected in a recent highly-cited paper on oscillatory flows of 
second-grade fluids [Fetecau & Fetecau (2005). Int. J. Eng. Sci., 43, 781-789]. The corrected solutions are shown 
to agree identically with numerical ones generated by a finite-difference scheme, while the original ones of Fetecau 
& Fetecau do not. A list of other recent papers in the literature that commit the error corrected in this Comment is 
compiled. Finally, a summary of related erroneous papers in this journal is presented as an Appendix. 

Keywords: Stokes' second problem (transient form), Second-grade fluid, Integral transforms 



1. Introduction 



Fetecau & Fetec au (2005) considered the transient version of Stokes' second problem for a second-grade (SG) 
fluid. Unfortunately, due to a mathematical error in imposing the boundary condition, an incorrect and unphysical 
solution is obtained. The purpose of the present Comment is to correct this error, which has unfortunately been 
promulgated through scores of papers, and to offer a new perspective on the subject. 

In the literature, many errors (of increasing variety and frequency) have been made by a growing list of authors 
attempting to obtain solutions to initial-boundary-value problems arising from simple flows of non-Newtonian fluids. 
For i n-depth discussions (including identification of the mistakes and how to correc t t hem) in a variety of contexts see , 
e.g.. jjordan et all 12004 1 Jordan & : Puril|2004|Jordaill2005l|Panto]gato & Jordani 120091; [Rees , 

2009llJordanll201CtlChristov & Christovil201ot IChristovL 12 010; An dersson & KumaranLuOlOtlLipscombeLuOlOtlSohl 



2010[|Munawar et al.ll201 lHChristovll201 it iKoshelL |2012F_ 

Another unsuccessful attempt, concurrent with that oflFetecau & Fetecau d2005h . at solving the transient version 
of Stokes' second problem for a SG fluid can be found in (lAsghar et al.Ll2006l). We s hall not concern ourselves with the 
latter work in the discussion below because, it suffices to note. lAsghar et al. N2006I) solve a linear constant-coefficient 
second-order ordinary differential equation by a perturbation series for which they take a dimensional quantity as the 
small parameter without giving a point of reference. 
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2. Start-up problem for a half-space 



Fetecau & Fetecaul (120051) first attempt to solve the following initial-boundary-value problem (IB VP): 



du 



d 2 u d 3 u 
dy 2 dydtdy ' 



u(0,t) = U f(<ot)H(t), 
u(y, t) — » as y — » 
u(y, 0) = 0, 



(y, t) € (0, oo) x (0, oo), 

t > 0, 
t > 0, 
y > 0. 



(la) 

(lb) 
(lc) 
(Id) 



Here, we break with the notation of iFetecau & Fetecaul d2005l) because of the similarity of the symbols for viscosity 
and velocity employed therein. To this end, v is the kinematic viscosity, a is the ratio of the first material modulus of 
the SG fluid to its density, u = u(y, t) is the fluid velocity in the x-direction, Uo and a> are, respectively, the amplitude 
and the frequency of the oscillations of the infinite plate upon which this SG fluid rests. Two cases are considered: 
/(■) = cos(-)and/(-) = sin(-). 



The Heav iside unit step function H{t) multiplying Uof(a>t) in Eq. dTbl is neglected bv lFetecau & Fetecaul d2005l) . 
Stokesl dl85ll p. 101-102), though lacking the mathematical tools such as the Heaviside function or the operational 
calculus, was able to clearly (and elloquently) explain the difference between imposing m(0, t) = Uof(a>t)H(t) and 
imposing u(Q, f) = Uof(ojt) for a Newtoni an fluid. Hence, though it is claimed that the problem in which "At time 
t — + the rigid plate begins to oscillate." ( Fetecau & Fetecaul 2005, p. 7 83) is solved, this is not the case. In fact, 
it is unclear whether the expressions derived bv lFetecau & Fetecau (2005) solve any physically relevant problem; for 
Stokes' first problem for a SG fluid, the similarly incorrect solutio n could be re-interpreted a s one arising from a 
fictitious body force or a different, but physical, boundary condition ( Christov & Christov . 2010l) . 



Upon applying the Fourier sine transform to Eq. dTal ). the second term on the right-hand side becomes 

dt 1 S^' Sin( ^ } dy = dt l-? u ^ + ?"(0, 0] = + f ^(0, 0, (2) 

where the first equality follows from assuming that u is twice continuously differentiable with respect to y on (0, oo), 
and that u, du/dy and d 2 u/dy 2 are integrable functions of y on (0, oo). Note that the decay condition stated in Eq. dTcb 
is short-hand notation for this. Here, an "s" subscript denotes the image of a function in the Fourier sine transform 
d omain, and £ denotes the Fo urier sine transform parameter. We have employed the Fourier sine transform convention 
of lChurchill & Brownl dl978l) . 



2.1. Case I: /(■) = cos(-). 

Using the Schwartz-Sobolev theory of distributions, also known as generalized functions (see, e.g.. lKanwallll983 
Chap. 2), we have, from Eq. dTbl . that 



^(0, t) = U f(tot)6(t) + coU f'(a>t)H(t) = U cos(tot)6(t) - coU sin(cjt)H(t) 
ot 



(3) 



for this boundary conditio n, where <5(-) is the Dirac delta distribution. By comparing Eq. OJ to the right-hand side of 
dFetecau & Fetecaul 120051 Eq. (3.4)), it is apparent that the error committed by the latter authors boils down to taking 
6(t) = and H(t) = 1. However, this is both physically and mathematically wrong as d(-) and H(-) are distributions, 
hence they do not have point va lues, and they are certa inly not equal to and 1, respectively, for all t. [A similar 
incorrect claim may be found in dTan & Masuoka . l2009l Eqs. (3) and (4)).] The equality in Eq. (O is not understood 
in a pointwise sense for individual values of t, rather it is understood as being true upon multiplication of both sides 
by a smooth function and inte gration over t e (0. oo) 

Another way to see wh y dTan & Masuokal 20091 Eqs. (3) and (4)), which can be considered as the reasoning 
behind the mistake made bv lFetecau & Fetecaul (120051) . are incorrect is to perform the following thought experiment. 



Let us specify, and we are free to do so, that the plate is set in motion at some t = tt, where fo > 0, while being at rest 
for all t < to. Then, all terms involving the Heaviside and Dirac distributions that appear in this new IB VP (whose 
solution we shall denote by 5) are of the form H(t - to) and 6{t - fo) with t € (0, oo) as before. Since we require 
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?o > 0, it should now be even more obvious that the claim 5(t - to) — and H(t — to) = 1 for all t e (0, oo) is incorrect. 
Finally, from the change of variable / i-» t + tp, it follows tha t u(y, t) - u(y, t + to) satisfies the original IB VP (Q~|i, 
demon strating that the claims made by iTan & Masuokal J2009h (and equivalent manipulations by IFetecau & Fetecaul 
(120051) ^ are erroneous. 

Completing the application of the Fourier sine transform in y and then applying the Laplace transform in t, the 
solution to the subsidiary equation in the Fourier-Laplace domain is 



k s (£ s) 



Uot 



s(l + ff£ 2 ) + v£ 2 



+ a 



(4) 



Note that this expression (correctly) takes into account all information given in Eq. (Q}. Here, a bar over a quantity 
denotes its image in the Laplace transform domain and s denotes the Laplace transform parameter. 

Inverting Eq. (0]i using partial fractions, a standard table of Laplace inverses and the definition of the Fourier sine 
transform, we arrive at 



u(y,t) = UoH(t) 



2 r 

* Jo 



£sin(;y£) 



v 2 £ 4 + w 2 (l + a^f 



v 2 e 
i + 



exp 



+ vcu sin(wf) + [v 2 f + au?(\ + afj\ cos(wf)| d^j . (5) 



We leave it to the reader to find the correct long-time (post-transient) expression for the solution. Here it is important 
to note that this limi ting form of the so l ution is not independent of time (i.e., it is not, strictly speaking, a "steady 
state," as claimed in (IFetecau & Fetecaul 120051) V 



2.2. Case II: /(■) = sin(-). 
Now, Eq. (0 becomes 



dt 



(0, t) = U sm(o)t)5(t) + coUo cos(cot)H(t). 



(6) 



In this case, however, there is a "self-canceling error," as discussed in dChristov & Christovll2010l:IChristovl l2010l). be- 
cause _£{sin(wf)<5(0) = 0. In other words, the term in the subsidiary equation that Fetecau & Fetecaul ( 2005I) neglected 
does not contribute to the solution. 

For this boundary condition, the solution of the subsidiary equation in the Fourier-Laplace domain is 



«s(£, S) = 



U { 



s(l + of 2 ) + v£ : 



lv<jj + aa>s \ 
s 2 + cj 2 )' 



(7) 



Inverting Eq. © using partial fractions, a standard table of Laplace inverses and the definition of the Fourier sine 
transform, we arrive at 



u(y,t) = UoH(t)\ 



2 -r- 2 

x Jo v- 



£sin(y£) 



e+com+ae) 2 \ v<JJSxp \ \+af 



>cos(cot) + [v^f + aa/(l + af)] sin(wf)| dA (8) 



Desp ite the "self-canceling error," we give this solution for completeness because it differs from (IFetecau & Fetecau , 



2005, Eq. (3.9)), perhaps due to a choice made during the inversion procedure. 
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3. Start-up problem for a strip 



Next, we turn to the problem considered in (IFetecau & Fetecaul 120051, Sect. 4). Restricting Eq. CD to the strip 
< y < d, we obtain a new IB VP: 



du d 2 u d^u 
~di = V df + a dydtdy' 

u(0, t) = U f(tot)H(t), 

u(d, t) = 0, 

u(y, 0) = 0, 



(y, t) e (0, oo) x (0, d), 

t>0, 
t > 0, 
ye(Q,d). 



(9a) 

(9b) 
(9c) 
(9d) 



Rather than follow the approach in (IFetecau & Fetecaul 120051 Sect. 4), here we make use of a significan tly more 
elegant solution procedure involving only the Laplace transform in time only as suggested bylJordarifeoOSl) . 

An eigenfunction expansion for a start-up problem (as attempted by IFetecau & Fetecaul d2005l) and some of the 
papers discussed in the Appendix of this Comment) presents many opportunities for making a mistake. If said eigen- 
function expansion were performed correctly, it would produce an identical solution to the one obtained by the more 
robust method we now employ. To this end, applying only the Laplace transform in t to Eqs. (f9ab— (f9cb. using the 
initial condition in Eq. (l9dt , and then solving the resulting subsidiary equation, we obtain 



u(y, s) 



U sinh [(l - lyl y/a) yjs/(l 2 K + si 2 )] 
sinh [ yjs/(l 2 K + si 2 )] 



x < 



/(•) = cos(-) 
/(•) = sin(-) 



(10) 



where we have introduced the dimensionless parameter / := d~ l -\[a and defined k :- v/a for convenience. 

To make the expression for the solution more tractable, we introduce the dimensionless independent v ariables 
y - y / yfa and t — Kt together with the dimensionless frequency Q. := oj/k. Now, it is readily established from ( Jordanl 
20051 Eq. (3.4)) and the convolution theorem for the Laplace transform, that the exact solution of IB VP © for the 
case /(■) = sin(-) is given by 



u(y, t) = U H(t) I sin(wf) (1 - ly) + ^ J] 



<t„ exp(-<T„f) - <x„ cos(Qf) - D. sin(Qf) 



(t 2 +Ql 2 



sin{lnny) 
n{\ + l 2 n 2 n 2 ) \ ' 



(11) 



where cr„ := l 2 n 2 n 2 (l +l 2 n 2 n 2 ) l . In turn, using Eq. (fTTT). the fact that Jj l {s) — <5'(t), and, once again, the convolution 
theorem for the Laplace transform, the exact solution of IB VP (O for the case /(•) = cos(-) is found to be 



( 2 °° 

u(y, t) = U H(t) | cos(wr) (1 - ly) - - J] 



cr 2 exp(-cr,,f) - cr„Q sin(Qf) + Q 2 cos(Of) 



oi + Q 2 



sin(ln7Ty) 
n(l+l 2 n 2 n 2 ) 



(12) 



Finally, we refer the reader to I Jordanl ( 2005 ) for a discussion of the equivalent error for the case of Stokes' first 
problem (see also Items[T]and|2]in the Appendix). 



4. Illustrating examples 



Us ing the unconditionally stable and formally second-order accurate finite-difference scheme o flChristov & Christov 
d2010l) . w e solve the IB VPs CD and Q numerically and compare the resulting solutions to the wrong transform so- 
lutions of IFetecau & Fetecaul (12005) and the corrected transfo rm solutions presented in Sects. [2] and [3] above. The 
scheme parameters are as given in Jchristov & Christov , 20ld Sect. 3). The integral representations of the analytical 
solutions are evaluated using the high-precision numerical integration routine Nlntegrate of the software package 
Mathematica (ver. 7.0.1). 

FigureQ~]shows a comparison of the erroneous, corrected and numerical solutions for the half-space problem with 
/(■) = cos(-). Similarly, Fig. [2] shows a comparison for the problem on a strip with /(■) = cos(-). We do not give 
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plots for the case of /(■) = sin(-) because the aforementioned self-canceling error renders the solution curves visually 
indistinguishable. 



Here, it is important to note another problem with the work of lFetecau & Fetecaul (|2005): when presenting numer- 
ical values of physical quantities units are neglected; a dditionally, values of quest ionable relevance are chosen for the 
material parameters. For example, for the purposes of (IFetecau & Fetecaul |2005[ Fig. 1), the viscosity and density of 
the fluid are taken to be those of glycerin, which (under normal conditions) cannot be claimed to be a SG fluid without 
justification. Arguably, it is impossible to pick a realistic value for a because undisputed experimental measurements 
do not exist ( Dunn & Raiagopaj. 1995 ). However, for the purposes of making a plot of the solution, w e can scale both 
v and a out by using the dimensionless variables y = y/ s/a, t = t{v/a), Q. = a)(a/v) and u - u/Uo (Bandelli et al 



1995). 



From Fig.[T] it is clear that the solution from Eq. d5 i agre es identically with the corresponding numerical solution. 

20051 Eqs. (3.6)) does not even satisfy th e boundary condition 



Meanwhile, the incorrect solution (IFetecau & Fetecau , 

rim v -»a+ u(y, t) = E/n (t > 0). Unlike the case of Stokes' first problem for a SG fluid corrected bv lChristov & Christov 



d2010h . ] lere it is unclear what kind of boundary condition the wro ng solution satisfies , or wh ether it has any physical 
meaning. Also, note the erroneous prediction, by the solution of lFetecau & Fetecaul (12.0051) . of a back-flow for y > 
2.25, preventing the solution from satisfying the asymptotic boundary condition in Eq. (fTctFI What is more, for 
the case of a strip, the curve drawn based on dFetecau & Fetecau , 2005 , Eq. (4.6)) appears to satisfy the boundary 
condition but exhibits highly unphysical oscillations near the moving plate. This is similar to the behavior of the 
incorrect solution to Stokes' first problem on a strip corrected bv lJorda nf (l2005b. 



For all cases considered in this section, the erroneous solutions from (IFetecau & Fetecaul 120 05) appear to even- 
tually agree with the correct an d numerical solutions in the limit of i » 1, though we are unable to confirm the 
long-time expressions giv en in (IFetecau & Fetecaul 120051) . In addition, for br evity, we did not check th e results in 
(IFetecau & Fetecaul |2005[ Sect. 5) for correctness. Moreover, any equation in (IFetecau & Fetecaul 120051) that is not 
explicitly corrected in the discussion above should not be assumed to be correct. 




(a)r = 1 




(b)t= 3.5 



Figure 1: (Color online.) u vs. y for the half -plane problem with /(•) = cos(-) and Q = 2n. Dashed: incorrect 
solution ( Fetecau & Fetecaul 2005 . Eq. (3 .6)); solid: correct solution [Eq. (0 above]; dots: numerical solution [using 
the scheme in dChristov & Christovll2010l Eq. (4))]. 



5. Conclusion 

The contribution of the present work is in correcting a number of mistakes made in an attempt to solve the tran- 
sient version of Stokes' second problem for a SG fl uid. Where possible , we ha ve also suggested alternative analytical 
approaches, which prevent the errors committed by IFetecau & Fetecaul (120051) . Unfortunately, the very same mistake 



: 'Pantokratoras (2008a b) uncovered a number of papers on steady boundary layer flows of non-Newtonian fluids in which the supposed solution 
also exhibits this kind of obviously unphysical behavior. 
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Figure 2: (Color online.) u vs. y for the strip problem with 1=1, /(•) = cos(-) and Q = 2n. Dashed: 100 terms of 
the incorrect series solution dFetecau & Fetecaul 12005 , Eq. (4.6)); solid: 100 term s of the correct solution [Eq. ( TT2l 
above]; dots: numerical solution [using the scheme in (Ichristov & ChristovLl20ld Eq. (4))]. 



has al re ady been committed i n the following p a pers: lAksel et alJ (120061): Fetecau & Fetecau (2006a); Kha n et al 
(|2007[) : iFetecau et all d2008h : iKhan et all J2008i l2009h : iHussain et all d2010h : iKhan et all d2010blah : lAhmad et al 
(l2010h : lYao&Liul(l2010l) . 

The present Comment is also useful when faced with a supposed solution to unidirectional flow of a SG, Oldroyd- 
B or the so-called "generalized" SG and Oldroyd-B fluids executing the same motion. Specifically, there is a parameter 
regime in which the SG flui d, studied herein under the transient version of Stokes' second problem, and studied in 
( Christov & Christov . 2010h . under Stokes' first problem, is reproduced. Hence, authors of such studies must be 
abl e to show their Fourier-Lap lace domain solution agrees (in the proper parameter limit) with, e.g., Eq. (0), Eq. Q 
or ( Christov & Christov , 2010l Eq. (2)). It is worth noting, ho wever, that not a ll derivations leading to the latter 
Fourier-Laplace domain solutions are mathematically correct (see lChristovil201 lb . 



Acknowledgments 

The authors would like to thank Prof. C. I. Christov for his advice and encouragement, and Prof. K. R. Rajagopal 
for informative discussions. The time-shift argument presented in Sect. I2.1l was suggested by one of the anonymous 
reviewers. 



References 

Ahmad, M., Zaman, H., & Rehman, N. (2010). Effects of hall current on unsteady MHD flows of a second grade fluid. Cent. Eur. J. Phys., 8, 
422^431. 

Aksel, N., Fetecau, C, & Scholle, M. (2006). Starting solutions for some unsteady unidirectional flows of Oldroyd-B fluids. Z. angew. Math. Phys. 
(ZAMP),57, 815-831. 

Andersson, H. I., & Kumaran, V. (2010). Comment on "Unsteady flow of a second grade fluid film over an unsteady stretching sheet" [Math. 

Comput. Modelling 48 (2008) 518-526]. Math. Comput. Modelling, 52, 1706-1707. 
Asghar, S., Nadeem, S., Hanif, K., & Hayat, T. (2006). Analytic solution of Stokes second problem for second-grade fluid. Math. Proh. Eng. ,2006, 

72468. 

Bandelli, R., Rajagopal, K. R., & Galdi, G. R (1995). On some unsteady motions of fluids of second grade. Arch. Mech., 47, 661-667. 
Christov, C. I., & Jordan, R M. (2009). Comment on "Stokes' first problem for an Oldroyd-B fluid in a porous half space" [Phys. Fluids 17, 023101 
(2005)]. Phys. Fluids, 21, 069101. 

Christov, I. C. (2010). Stokes' first problem for some non-Newtonian fluids: Results and mistakes. Mech. Res. Commun., 37, 717-723. 
larXiv: 1009.44161 

Christov, I. C. (2011). Comments on: "Energetic balance for the Rayleigh-Stokes problem of an Oldroyd-B fluid" [Nonlinear Anal. RWA 12 

(2011) 1]. Nonlinear Anal. RWA, 12, 3687-3690. larXi v : 1 107 . 2547 1 
Christov, I. C, & Christov, C. I. (2010). Comment on "On a class of exact solutions of the equations of motion of a second grade fluid" by C. 

Fetecau and J. Zierep (Acta Mech. 150, 135-138, 2001). Acta Mech., 215, 25-28. larXiv: 1003.21881 
Churchill, R. V., & Brown, J. W. (1978). Fourier Series and Boundary Value Problems. New York: McGraw-Hill. 



6 



Dunn, J. E., & Rajagopal, K. R. (1995). Fluids of differential type: Critical review and thermodynamic analysis. Int. J. Eng. Sci., 33, 689-729. 

Erdogan, M. E., & Imrak, C. E. (2007). On the comparison of the methods used for the solutions of the governing equation for unsteady unidirec- 
tional flows of second grade fluids. Int. J. Eng. Sci., 45, 786-796. 

Fetecau, C, & Fetecau, C. (2005). Starting solutions for some unsteady unidirectional flows of a second grade fluid. Int. J. Eng. Sci., 43, 781-789. 

Fetecau, C, & Fetecau, C. (2006a). Stalling solutions for some simple oscillating motions of second-grade fluids. Math. Prob. Eng., 2006, 23587. 

Fetecau, C, & Fetecau, C. (2006b). Starting solutions for the motion of a second grade fluid due to longitudinal and torsional oscillations of a 
circular cylinder. Int. J. Eng. Sci., 44, 788-796. 

Fetecau, C, Hayat, T., Fetecau, C, & Ali, N. (2008). Unsteady flow of a second grade fluid between two side walls perpendicular to a plate. 
Nonlinear Anal. RWA, 9, 1236-1252. 

Hayat, T., Asghar, S., & Siddiqui, A. (2000). Some unsteady unidirectional flows of a non-Newtonian fluid. Int. J. Eng. Sci., 38, 337-346. 

Hayat, T., Siddiqui, A., & Asghar, S. (2001). Some simple flows of an Oldroyd-B fluid. Int. J. Eng. Sci., 39, 135-147. 

Hayat, T., Siddiqui, A., & Asghar, S. (2005). Exact solutions for the unsteady rotational flow of non-Newtonian fluid in an annular pipe. Int. J. 
Eng. Sci., 43, 281-289. 

Hussain, M., Hayat, T., Asghar, S., & Fetecau, C. (2010). Oscillatory flows of second grade fluid in a porous space. Nonlinear Anal. RWA, 11, 
2403-2414. 

Jordan, P. M. (2005). A note on start-up, plane Couette flow involving second-grade fluids. Math. Prob. Eng., 2005, 539-545. 

Jordan, P. M. (2010). Comments on: "Exact solution of Stokes' first problem for heated generalized Burgers' fluid in a porous half-space" 

[Nonlinear Anal. RWA 9 (2008) 1628]. Nonlinear Anal. RWA, 11, 1198-1200. 
Jordan, P. M., & Puri, A. (2004). Noncausal effects arising from a Maxwell fluid model with negative relaxation. Eur. J. Phys., 25, 829-834. 
Jordan, P. M., Puri, A., & Boros, G. (2004). On a new exact solution to Stokes' first problem for Maxwell fluids. Int. J. Non-Linear Mech., 39, 

1371-1377. 

Kanwal, R. P. (1983). Generalized Functions: Theory and Technique. New York: Academic Press. 

Khan, M., Ali, S. H., & Qi, H. (2009). Exact solutions for some oscillating flows of a second grade fluid with a fractional derivative model. Math. 
Comput. Model, 49, 1519-1530. 

Khan, M., Anjum, A., Fetecau, C, & Qi, H. (2010a). Exact solutions for some oscillating motions of a fractional Burgers' fluid. Math. Comput. 
Model, 51, 682-692. 

Khan, M., Anjum, A., Qi, H., & Fetecau, C. (2010b). On exact solutions for some oscillating motions of a generalized Oldroyd-B fluid. Z angew. 
Math. Phys. {TAMP), 61, 133-145. 

Khan, M., Naheed, E., Fetecau, C, & Hayat, T. (2008). Exact solutions of starting flows for second grade fluid in a porous medium. Int. J. 
Non-Linear Mech., 43, 868-879. 

Khan, M., Saleem, M., Fetecau, C, & Hayat, T. (2007). Transient oscillatory and constantly accelerated non-Newtonian flow in a porous medium. 

Int. ]. Non-linear Mech., 42, 1224-1239. 
Koshel, K. (2012). Comments on "Peristaltic flow of a Williamson fluid in an asymmetric channel" (Commun. Nonlinear Sci. Numer. Simul. 15, 

1705-1716). Commun. Nonlinear Sci. Numer. Simul, 17, 483^184. 
Lipscombe, T. C. (2010). Comment on 'Application of the homotopy method for analytical solution of non-Newtonian channel flows'. Phys. Scr., 

81, 037001. 

Munawar, S., Mehmood, A., & Ali, A. (201 1). Comment on "Analysing flow and heat transfer of a viscoelastic fluid over a semi-infinite horizontal 

moving flat plate, IJNLM, 43 (2008) 772". Int. J. Non-Linear Mech., 46, 1280-1282. 
Pantokratoras, A. (2008a). Comment on "Similarity analysis in magnetohydrodynamics: Hall effects on forced convective heat and mass transfer 

of non-Newtonian power law fluids past a semi-infinite vertical flat plate" by A. A. Afify, E.M. Aboeldahab and E.S. Mohamed (Acta Mech. 

177, 71-87, 2005). Acta Mech., 198, 119-121. 
Pantokratoras, A. (2008b). Comment on "Thermal radiation effects on MHD flow of a micropolar fluid over a stretching surface with variable 

thermal conductivity", by Mostafa A.A. Mahmoud [Physica A 375 (2007) 401-410]. Physica A, 387, 2416-2418. 
Rees, D. A. S. (2009). Comment on "Effects of chemical reaction on free convective flow of a polar fluid through a porous medium in the presence 

of internal heat generation" by Patil and Kulkarni, Int. J. Thermal Sci. 47 (2008) 1043-1054. Int. J. Thermal Sci., 48, 847-848. 
Soh, C. W. (2010). Comment on 'A point source solution for unidirectional flow of a viscoelastic fluid' [Phys. Lett. A 372 (2008) 4041]. Phys. 

Lett. A, 374, 2098-2100. 

Stokes, G. G. (1851). On the effect of the internal friction of fluids on the motion of pendulums. Trans. Cambridge Phil. Soc, 9 (Part II), 8-106. 

http : //www .biodiversitylibrary . org/item/19878#page/216/ 
Tan, W., & Masuoka, T. (2009). Response to "Comment on 'Stokes' first problem for an Oldroyd-B fluid in a porous half space'" [Phys. Fluids 21, 

069101 (2009)]. Phys. Fluids, 21, 069102. 
Yao, Y, & Liu, Y. (2010). Some unsteady flows of a second grade fluid over a plane wall. Nonlinear Anal. RWA, 11, 4442^1450. 
Zierep, J., & Fetecau, C. (2007). Energetic balance for the Rayleigh-Stokes problem of a second grade fluid. Int. J. Eng. Sci., 45, 155-162. 



Appendix. Other erroneous recent papers on simple flows of non-Newtonian fluids in this journal 

1. lHavat et al.l (l2000l) present, in their Sect. 2.3, an incorrec t solution for th e plane Cou e tte flo w of a SG fluid 
between parallel plates. The correct solution is given by Jordan ( 2005 ). Hayat et al. ( 2000l) have neglected 
to include H{f) in their initial condition in their Eq. (12), which renders the transformation in their Eq. (13) 
erroneous. If the H(t) pre-factor were carried through, there would be an additional term proportional to S(t) 
on the left-hand side of their Eq. (15 ) . Wh at is more, the initial condition, which must be understood as the 
limit t — > 0~ of the solution dJordanl EoToh . would require that the right-hand side of their Eq. (18) is zero, 
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identically, even after the attempted transformation. The error in ( IHavat et al. . 2000. Sect. 2.3) also impacts a 
number of results given in their Section 3. 

Havat et al.l (120011) present, in their Sect. 6, incorrect solutions for the plane Couet te flow of an Oldroyd -B fluid 



between parallel plates. The mistake is the same as described in Item[T] again, see Jordan (2005, 2010h . 
Fetecau & Fetecau d2006bl) omit H(t) in their Eq. (7). As in Item Q] above, this leads to the lack of a term 
proportional to 6{t) on the left-hand side of their Eq. (10.1), and the initial condition in their Eq. (10.3) is 
incorrect. Therefore, when the transformation given in their Eq. (9) is made, it is implicitly assumed the outer 
cylinder has been and will continue to oscillate. Hence, the problem solved is not that of start-up as claimed. 
Erdogan & Imrakl ( 2007 ) present, in their Sect. 4, incorrect solutions to and incorrect discussion of Stokes' first 



problem for a SG fluid. These are corrected in dChristov & ChristoviuOlOl) . 



Havat et al.l ( 12005b claim that Sect. 3 of their paper presents "Poiseuille flow due to a constant pressure gradient." 



However, the start-up flow for a cylinder that begins to rotate at constant angular velocity is considered instead. 
Unfortunately, H(t) is missing from the boundary condition in their Eq. (6), which makes the transformation in 
their Eq. (7) and the resulting governing Eq. (8) incorrect. This is the same mistake as in Item[3]above. 



Zierep & Fetec au (2007) perform an energy analy sis using an incorrect soluti on, rendering their analysis erro- 



neous. The correct solution to be used is given in dChristov & Christovll20ld Eq. (3)) 
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